Effect of a magnetic flux on the critical behavior of a system with long range hopping 
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We study the effect of a magnetic flux in a ID disordered wire with long range hopping. It is shown 
that this model is at the metal-insulator transition (MIT) for all disorder values and the spectral 
correlations are given by critical statistics. In the weak disorder regime a smooth transition between 
orthogonal and unitary symmetry is observed as the flux strength increases. By contrast, in the 
strong disorder regime the spectral correlations are almost flux independent. It is also conjectured 
that the two level correlation function for arbitrary flux is given by the dynamical density-density 
correlations of the Calogero-Sutherland (CS) model at finite temperature. Finally we describe the 
classical dynamics of the model and its relevance to quantum chaos. 
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The addition of disorder to an otherwise metallic sam- 
ple strongly modifies its properties. As disorder increases 
eigenf unctions start to localize. A MIT is observed in 
the thermodynamic limit in systems of dimension greater 
than two. The moments of the eigenfunctions at the 
MIT show anomalous scaling with the sample size [1] L, 
P q = J d d r\ip(r)\ 2q oc L~ D i( q ~ 1 \ P q is the inverse par- 
ticipation ratio and D q is a set of exponents describing 
the transition. Thus the scaling at the MIT is in between 
that of a good metal, P q oc L - "^' -1 ), and that of an insu- 
lator (localized eigenfunctions), P q oc L°. Eigenfunctions 
with such anomalous scaling are named "multifractals" , 
for a review see [2]. Signatures of a MIT are found not 
only in the eigenfunctions [3] but also in the spectral 
fluctuations [4]. It was argued in [5] that although level 
repulsion typical of the metallic regime is still present 
at the MIT, the long range correlations are weaken due 
to the multifractal character of the wavefunctions. The 
number variance was claimed to be asymptotically linear 
with a slope around 0.25 [5]. These features, level re- 
pulsion and sub-Poisson number variance combined with 
the scale invariance of the spectrum [4] are named 'criti- 
cal statistics' [6] and are considered genuine spectral sig- 
natures of a MIT. Eigenfunction and spectral properties 
are indeed related, for multifractal eigenstates not too 
sparse the slope of the number variance x is related to 
the multifractal exponent Di by \ = d ~2d 2 \J\ > wnere d 1S 
the dimension of the space. Different generalized random 
matrix model (gRMM) have been successfully employed 
to describe 'critical statistics' [8]. Recently [9,10], the 
spectral correlations of one of those gRMM was demon- 
strated to be equivalent to the spatial correlations of the 
Calogero-Sutherland (CS) model [11] at finite tempera- 
ture. Temperature in the CS model is related with the 
multifractal exponent D 2 at the transition. 

Short range Anderson [12] models have been broadly 
utilized for numerical investigation at MIT. However in 
certain systems, like glasses with strong dipole interac- 
tions [15], long range hopping is possible and thus the 
Anderson model must be modified accordingly. Although 
the introduction of long range hopping dates back to the 



famous Anderson's paper on localization [12], these mod- 
els did not attract too much attention until the numeri- 
cal work of Oono [13,14] and the renormalization group 
treatment of Levitov [15] in the context of glassy systems. 
The main conclusion of these works was that power low 
hopping may induce a MIT even in one dimensional sys- 
tems if the exponent of the hopping decay matched the 
dimension of the space. A related problem, a random 
banded matrix with a power law decay was discussed 
in [16]. It was shown analytically that for a 1/r band 
decay the eigenstates are multifractal and the spectral 
correlations resemble those at the MIT. Intense numeri- 
cal study in recent years [17] has corroborated the close 
relation between this random banded matrix model and 
the Anderson model at the MIT. 

Another related issue of current interest is the effect 
of a magnetic field at the MIT [18,19]. In the metal- 
lic regime, an exact analytical treatment was developed 
in [20]. It turns out that, in this limit, the two level 
correlation function describing the crossover between or- 
thogonal (no flux) and unitary symmetry is equivalent 
to the dynamical density-density correlations of the CS 
model [21] (see also [22] for different initial conditions) 
provided that "time" in the CS model is traded by mag- 
netic field in the disordered system. Numerical results at 
the 3D MIT [18] show that the magnetic flux still influ- 
ences the short-range spectral correlations. The impact 
on long range correlations is still not settled [19] though 
there is agreement that the effect, if any, must be small. 

One of the main aims of this paper is to further clarify 
this issue by investigating a one dimensional (ID) wire 
with power low hopping and a magnetic flux attached to 
it. This model presents typical features of a MIT for all 
values of disorder with the advantage that much larger 
volumes can be simulated. We also conjecture, based on 
the analogy with the CS model above mentioned, a exact 
expression for the two level spectral correlation function 
at the MIT for arbitrary flux. 

Our starting point is the following Hamiltonian de- 
scribing a ID closed wire with long-range hopping, 
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where, n, m label the lattice sites, r = n — m, e n and e„ m 
are two sets of random number both distributed in a box 
[-W, W], F(r, b) is the long range hopping term, 
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< A < oo is the band size and N is the system size Since 
F(r, b) ~ 1/r for N r our model is critical, namely re- 
produce typical features of a MIT for all values of A [16]. 
Sine-like interaction is introduced in order to assure pe- 
riodicity. The phase factor (Peierls substitution) above 
is related to the magnetic flux (in units of the funda- 
mental flux h/c) by, 8 nm = f Adr = a(n — m)/N for 
n — m< N/2 and 8 nm = J Adr = a(l — (n — m) /N) for 

n — m> N/2, a is the magnetic flux. A is the (constant) 
potential vector. Such Hamiltonian may be relevant for 
glassy systems [14], quantum effect of classical anomalous 
diffusion [23] and quantum chaos [25]. 

The spectral fluctuations are studied by direct diago- 
nalization of the Hamiltonian (1) for different sizes rang- 
ing from N = 500 to N = 8000. The eigenvalues thus 
obtained are unfolded with respect to the mean spectral 
density. The number of different realizations of disor- 
der is chosen such that for each N the total number of 
eigenvalue be at least 2 x 10 5 . In order to avoid edge 
effects, only 30% of the eigenvalues around the center of 
the band are utilized. Without loss of generality we have 
set W = 1. 

First we investigate short range correlations by eval- 
uating P(s), the level spacing distribution (LSD). This 
correlator gives the probability of having two eigenvalues 
at a distance s. In the insulator regime (uncorrelated 
eigenvalues) P(s) = e~ s , in the metallic regime (Wigner- 
Dyson statistics) P(s) ~ s@ for s <§C 1 (the presence of a 
flux drives the spectrum from = 1 orthogonal to /3 = 2 
unitary symmetry). At criticality one expects level repul- 
sion as in Wigner-Dyson statistics s«l and exponential 
decay for s ^> 1, P(s) ~ e~ KS with k > 1 [26]. In figure 1 
we plot P(s) for different flux values a and A = 2. As the 
flux strength increases, a transition is clearly observed be- 
tween orthogonal and unitary symmetry. For s <C 1 (see 
inset) level repulsion is still present, P(s) ~ s' 3 is Wigner- 
Dyson like but with the prefactor modified by the flux. 
The strong disorder regime is reached by choosing a cou- 
pling constant (A = 0.44) such that the resulting spec- 
trum resembles the one at the 3D MIT. As observed in 
figure 2, P(s) is almost independent of the flux strength. 
For s > 1 the decay is exponential (see inset) and inde- 
pendent of the flux. For s -C 1 , the effect of the flux is 
still important, only for a ~ the expected P(s) ~ s be- 
havior for orthogonal symmetry is recovered. Such result 
agrees with previous [18,19] numerical simulations. Let 
us now move to long range spectral correlations. 




FIG. 1. Level spacing distribution P(s) for X — 2 and dif- 
ferent fluxes a. A transition between the Gaussian Orthogonal 
(GOE) and Gaussian Unitary (GUE) ensemble is observed as 
the flux is increased, even the asymptotic behavior of P(s) is 
modified by the flux. 



The number variance S 2 (L) = (L) 2 — {L 2 ) = L — 
2 J Q ds(L — s)R2(s,a) (i?2(s,a) is the two level corre- 
lation function) measures the stiffness of the spectrum. 
Fluctuations are small in the metallic regime (Wigner- 
Dyson statistics) with £ 2 (£) ~ log(i) for L ^> 1. For 
localized eigenstates, the eigenvalues are uncorrelated 
(Poisson statistics) and S 2 (£) = L. As mentioned ear- 
lier, the number variance is asymptotically proportional 
to xL (x < 1) at the MIT. In figure 3 we plot the number 
variance in the weak disorder regime (A = 2) for different 
flux values. As the flux is increased, a smooth transition 
is observed between critical GOE [10,9] and critical GUE 
[8]. We remark the slope is not modified by the flux. 
This is in apparent disagreement with the results for a 
2D Anderson model in the weak multifractality regime 
where the slope is two times smaller in the unitary case 
[20]. The reason for that is the way in which the flux is 
introduced due to the long range hopping. If the flux in 
(1) is introduced as usual in short range models the 2D 
Anderson model results are recovered. The asymptoti- 
cally linear number variance (see fig 3) together with the 
level repulsion indicates that our model is described by 
critical statistics [6]. The strong disorder regime is ex- 
plored by choosing a coupling constant (A = 0.44) such 
that the slope of the number variance x ~ 0.29 be similar 
to that at the 3D MIT. We observe (fig 4) that the flux 
only modifies the spectral correlations up to (L ~ 1) (see 
inset fig 4). Beyond this point the number variance is 
independent of the flux, thus suggesting that symmetry 
is washed out by disorder. We want to stress that the 
transition from weak to strong disorder is not sharp. As 
A decreases, the flux impact decreases smoothly. Curi- 
ously, until values of A very close to those at the 3D MIT, 
the effect, although small, is still sizable. Such behavior 
makes more difficult an accurate account of the flux at 
the 3D MIT. Finally let us mention that although not 
shown in the figures it was explicitly checked that the 
spectral correlators are not dependent on the size of the 
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matrix as expected at a MIT. 




FIG. 2. P(s) for A = 0.44 and different fluxes. Only for 
s <C 1 the system is sensitive to the flux. For s 3> 1 all curves 
become indistinguishable (inset). 

After discussing the numerical results we propose 
an exact analytical relation for the two level cor- 
relation function of (1). We claim that such 
correlation function is equivalent to the dynamical 
density-density correlations of the CS model [11], 
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at finite temperature where the position of the CS par- 
ticles corresponds with the (unfolded) eigenenergies of 
(1), "time" ~ a 2 and temperature h is related to A. 
Let us first argue how these density-density correlations 
are obtained [10,21]. The density of probability of the 
CS ground state (zero temperature) for (/? = 1,2,4) 
is equivalent to the joint distribution of eigenvalue of 
the Gaussian random matrix models (GRMM). In the 
large N limit, the two point spectral correlations of the 
GRMM are explicitly expressed through a spectral ker- 
nel. Such kernel in the language of the CS model cor- 
responds with the amplitude of probability of having 
two free fermions in a ensemble of N around two arbi- 
trary positions. The idea is that, based on the Luttinger 
liquid nature of the CS model [24], the density-density 
correlations of the CS model at finite low temperature 
can be obtained by using the results of GRMM but re- 
placing the spectral kernel above mentioned by its finite 
temperature analogue. In the grand canonical ensemble 

[8], K(s = x- y,0) h = £„ g^gfel where V™ are 
the single particle wave functions for free fermions and 
E n = n/N. In the large N limit the wave functions are 
replaced by plane waves with energy given by k 2 . 

Dynamical density-density correlations [21] can be in- 
cluded in this formalism by expressing the time depen- 
dent density as p(s,t ~ a 2 ) = e Hcst p(s)e- Hcst . The 
crossover between orthogonal and unitary symmetry in 
the GRMM [20] corresponds in the language of the CS 
model with density-density correlations with initial con- 
ditions given by particles distributed according to the CS 
model for (3=1 and then evolving for t > according 
to the CS model for (3 = 2 (free fermions). The density- 



density correlations in this case [21] are also expressed 
through a spectral kernel as the one mentioned above but 
now involving the propagation of free fermions (hole and 
particle). Combining both effect, the dynamical (with 
the inital conditions above mentioned) density-density 
correlations of the CS model at finite temperature are 
given by (see [9,21] for details), 

(p(s, a)p(0, 0)) h - (p(s, a)) h (p(0, 0)) h = R 2 (s, a) h = (3) 
K 2 P (s, Q) h - (JfcKpfa a)h) J K H {s', a) h ds' 

r-OO 

I 2 2 
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where n#(fc) 
pation number, 
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is the hole occu- 
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z = l/(e 1 ^ h — 1) is the fugacity, a is the 
magnetic flux (time), and h is related to A by 2nh ~ 1/A 
for h< 1 [16]. 

The a — > limit in (3) yields the density-density corre- 
lations of the CS model at finite temperature h <C 1 [10,9] 
which resembles the two point spectral correlations of a 
disordered system at the localization-delocalization tran- 
sition. The limit h — > corresponds with the exact an- 
alytical result for the crossover between orthogonal and 
unitary symmetry of a disordered system in the metallic 
regime [20] which is equivalent to the exact dynamical 
density-density correlations of the CS model [21]. 
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FIG. 3. Number variance S 2 (L) for weak disorder A = 2 
and different flux strengths. Lines are the analytical prediction 
(3) and symbols are the results obtained from the numerical 
diagonalization of (1) for iV = 8000. 

The 8 > C 1 limit of (3) also agrees with the 
density-density correlations obtained in [24,22] by using 
conformal field techniques. 

Numerical calculations also support the validity of our 
conjecture. As observed in fig 3, the agreement between 
the analytical number variance based on (3) and the nu- 
merical result is excellent for the whole range of fluxes at 
A ~ 2. We remark that no fitting is involved as the value 
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of h utilized corresponds with the analytical prediction 
h = ~ 0.079 [16]. In the strong disorder regime (see 
figure 4) the agreement is also remarkable but in this 
case no relation between h in A is known so the param- 
eter h — 0.267 is the best fit to the numerical results. 
We want to point out that further work is needed to test 
whether the conjecture (3) is valid beyond the h <C 1 
limit. 



E (L) 




FIG. 4. Number variance E 2 (L) for a disorder (A = 0.44) 
mimicking a true 3D Anderson transition. The lines represent 
the analytical prediction (3) (h — 0.267) for no flux (segment) 
and maximum flux (solid) . The symbols correspond to the nu- 
merical diagonalization of (1) for N = 8000. The dependence 
on the flux observed for small L ~ 1 (see inset) is also repro- 
duced by (3). 

Finally we study the relation between classical and 
quantum properties of our model. For a = 0, h <C 1, 
the two level spectral correlation function (s 3> 1) can 
be expressed through the spectral determinant of a dif- 
fusion operator [28] describing the classical dynamics of 
(1), 

Ra(s) = - A^S^ + 2cos(2™)e aGf M (4) 

where, e G ^ = ,£)-!(«) = U m ^ is the s P ec " 

tral determinant and e m = m/h are the eigenvalues of 
the (anomalous) diffusion operator. As expected, the 
(s 3> l,a = 0), limit of (3) coincides with (4). Unlike 
short range models, the type of classical diffusion asso- 
ciated with (1) is anomalous [23]. Such diffusion is de- 
scribed by a fractional Fokker-Planck equation [29]. For 
band decaying as 1/r 7 , 1/2 < 7 < 3/2, the resulting clas- 
sical dynamics is superdiffusive with |r| ~ fV( 2 7-i) [16], 
In the our case (7 = 1) the associated fractional Fokker 
Planck equation is first order in space [30] and the eigen- 
values of the diffusion operator (4) are thus linear in m 
instead of quadratic as for normal diffusion. This link 
between anomalous classical motion and quantum corre- 
lations at the MIT may be utilized to find out conditions 
for the Anderson transition in quantum chaotic systems. 

In conclusion, we have performed a numerical and ana- 
lytical investigation in a ID wire with long range disorder 



and a flux attached to it. In the weak disorder regime 
we have observed that despite the critical character of 
the model the effect of the flux is still important. By 
contrast, in the strong disorder limit, the spectral cor- 
relations are almost independent of the flux except for 
small eigenvalue separations. We have also conjectured, 
by exploiting an analogy with the CS model, an exact an- 
alytical relation for the spectral correlations at the MIT 
valid for arbitrary flux and disorder We have sug- 

gested that, at least in ID, the MIT, a quantum mechanic 
phenomenon, may be related with certain features of the 
classical dynamics (anomalous diffusion). We expect this 
relation to be relevant in quantum chaos problems. 
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